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Evolution of solar-type stellar wind 
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By extending our self-consistent MHD simulations for the solar wind, we study the evolution of stellar winds of solar-type 
stars from early main sequence stage to red giant phase. Young solar- type stars are active and the mass loss rates are 
larger by up to ~ 100 times than that of the present-day sun. We investigate how the stellar wind is affected when the 
magnetic field strength and fluctuation amplitude at the photosphere increase. While the mass loss rate sensitively depends 
on the input energy from the surface because of the global instability related to the reflection and nonlinear dissipation 
of Alfvenwaves, it saturates at ~ 100 times because most of the energy is used up for the radiative losses rather than the 
kinetic energy of the wind. After the end of the main sequence phase when the stellar radius expands by ~ 10 times, the 
steady hot corona with temperature 10 6 K, suddenly disappears. Chromospheric materials, with hot bubbles embedded 
owing to thermal instability, directly stream out; the red giant wind is not a steady stream but structured outflow. 
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1 Introduction 

The origin of the energy which drives the solar wind is be- 
lieved to be in the surface convection zone. Turbulent mo- 
tion associated with the convection excites various modes of 
waves that propagate upwardly. Magnetic field lines, which 
are also a consequence of dynamo processes in the surface 
convective motion (e.g. Yoshimura 1975; Parker 1993), play 
as a guide for the upgoing waves. Among various types of 
waves, the Alfvenwave is the most reliable candidate which 
transfers the energy from the surface convection to the up- 
per region where the solar wind is accelerated because it can 
travel a longer distance due to the incompressive nature (e.g. 
Alazraki & Couturier 1971; Belcher 1971; Hollweg 1973). 

Alfvenwaves have been extensively studied in the con- 
text of driving the solar wind. The dissipation of Alfven 
waves is a key to understand the acceleration of the solar 
wind because it controls the energy and momentum trans- 
fer from the waves to the gas. Various types of dissipation 
mechanisms have been investigated such as ion-cyclotron 
resonance (e.g. Axford & McKenzie 1992; Kohl et al. 1998), 
turbulent cascade (e.g. Matthaeus et al.1996; Verdini & Velli 
2007), and the nonlinear mode conversion to compressive 
modes (e.g. Kudoh & Shibata 1999) by using steady-state 
models (e.g. Cranmer & van Ballegooijen 2007) or dynam- 
ical simulations (e.g. Suzuki & Inutsuka 2005). 

Among others, we have carried out magnetohydrody- 
namical (MHD) simulations covering from the photosphere 
to the sufficiently outer region where the solar wind is al- 
ready accelerated (Suzuki & Inutsuka 2005, 2006; Matsumoto 
& Suzuki 2012). One of the greatest advantage of our dy- 
namical simulations is that the inner boundary is not an ad 



hoc 'coronal base' or 'upper chromosphere' but the pho- 
tosphere where there are much observational information^ 
(e.g. Fujimura & Tsuneta 2009; Matsumoto & Kitai 2010). 
Therefore, mass loss rate can be directly determined as an 
output of the surface properties. 

This type of wind is not unique to the solar wind. Main 
sequence stars with mass comparable to or less than the 
solar mass possess surface convection zones, and the stel- 
lar winds emanating from these stars by similar mechanism 
(e.g. Cranmer & Saar 201 1). Observations of asterospheres 
of low-mass to intermediate-mass stars show that the mass 
loss rates of young stars are much larger up to ~ 100 times 
of the present solar value (Wood et al. 2002, 2005). 

Red giant stars also have surface convection zones. Alfven 
waves possibly play a role in driving the stellar winds from 
red giant stars (e.g. Hartmann & MacGregor 1980; Holzer 
et al.1983). 

In this paper, we discuss stellar winds originating from 
surface convection zones. In particular we focus on how the 
mass loss rate is determined from the photospheric proper- 
ties, based on our MHD simulations. 

2 Our previous works for the present-day 
solar wind 

In this section, we briefly summarize our MHD simulations 
for the solar wind. In Suzuki & Inutsuka (2005, 2006), we 
performed self-consistent one-dimensional (ID) simulations 
that handle the propagation, reflection, and dissipation of 
MHD waves in an individual super-radially open magnetic 
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1 Recent steady-state calculations (e.g. Verdini & Velli 2007; Cranmer 
& van Ballegooijen 2007) also include the photosphere. 
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flux tube from the photosphere to the solar wind region. We 
took into account radiative cooling and thermal conduction 
to examine the coronal heating without ad hoc assumptions 
and three components of magnetic field and velocity to treat 
Alfvenwaves. We injected transverse fluctuations from the 
photosphere and run the simulations until the quasi-steady 
states were achieved. A main result of these papers is that 
the input of the transverse fluctuations with ~ 1 km/s at 
the photosphere naturally drives the solar wind which is ob- 
served today. 

The Alfvenwaves which are excited from the photosphere 
are mostly reflected back downward because of the change 
of the Alfvenspeed (Moore 1991), but roughly 10 % of the 
initial Poynting flux associated with the Alfvenwaves pene- 
trates to the corona and contributes to the heating of the so- 
lar wind. The main channel of the dissipation of the Alfven 
waves is the generation of compressive waves, particularly 
slow MHD waves (Kudoh & Shibata 1999). The steepening 
of the wave fronts leads to the shock dissipation of these 
excited compressive waves. In general, however, in the ID 
simulations the shock dissipation tends to be overestimated 
because the waves are confined in the individual flux tubes 
without leakage. 

In Matstumoto & Suzuki (2012), we extended to two di- 
mensional (2D) simulations. The biggest difference is that 
we can treat the leakage of the waves to neighboring flux 
tubes. Cascading Alfvenicturbulence can be also handled, 
although it is restricted to the 2D space. In the 2D simula- 
tion, the effect of compressive waves is suppressed, while 
the turbulent cascade becomes comparably important in the 
dissipation of the Alfvenwaves. The overall dissipation rate 
of the Alfvenwaves is comparable to that of the ID cases, 
and as a result, the ID and 2D simulations give similar mass 
loss rates and terminal velocities. 

In the following sections, we extend these simulations 
from the present-day solar wind to more active solar-type 
stars (§3) and red giant stars (§4). Because the 2D and ID 
simulations yield similar global properties of the winds, we 
use the ID simulations to save computational time. 

3 Mass loss from young active solar-type 
stars 

Young solar-type stars are generally very active. The X-ray 
flux is up to ~ 1000 times larger than the present Sun (e.g. 
Giidel 2004), and the transition region flux is also much 
larger (Ayres 1997). Observation of young main sequence 
stars show very strong magnetic fields with an order of kG 
or even larger (e.g., Donati & Collier Cameron 1997; Saar & 
Brandenburg 1999). The mass loss rate is also much higher 
but seems saturated around ~ 100 times of the present solar 
level (Wood et al. 2002; 2005). 

In this section, we investigate how the solar atmosphere 
reacts to increases of magnetic field strength and fluctuation 
amplitude at the photoshere. As the reference case we set up 
a magnetic flux tube with field strength Bq = 1 kG at the 
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Fig. 1 Comparison of the wind structures. The solid lines 
are the wind structure for the active case (see text) and the 
dashed lines are the wind structure of the reference case. 
From top to bottom, velocity, density, and temperature are 
displayed. 

photosphere and a super-radial expansion factor / = 1000 
to match recent HINODE observations (Tsuneta et al. 2009; 
Ito et al. 2010), and inject fluctuations, Sv = 1.4 km s _1 
from the photosphere. This reference case well reproduces 
the average global properties of the present-day solar wind. 
We perform more than 30 simulation runs with different Bo, 
/, and Sv. 

Figure [TJ shows the stellar wind structure which adopts 
B = 2 kG, Sv = 2.8 km s" 1 , and / = 1000 (solid lines, 
labeled as 'Young Sun') in comparison with the reference 
case (dashed lines). This active case gives the mass loss 
rate, M = 4 x 10 -13 M©/yr, which is 20 times larger than 
2 x 10~ 14 M /yr obtained from the reference case (~ the 
present value). Since the injected Poynting flux (oc BqSv 2 ) 
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Fig. 2 Energetics of stellar wind derived from the simulations; each data point indicates each run. The left panel shows 
the fractions of the kinetic energy fluxes on the total input Poynting fluxes from the photospheres, and the right panel shows 
the fractions of the total radiative losses from the transition regions and coronae on the total input Poynting fluxes. 



in the active case is only 8 times larger, the mass loss rate 
quite sensitively depends on the energy input. This is be- 
cause of the global instability involved with the reflection 
and nonlinear dissipation of Alfvenwaves (see Suzuki 2012 
for details). 

The location of the transition region in the active case 
is at ~ l.li?©, which is much higher than the height of the 
transition region of the reference case (w 1.003i?©). This 
is because in the active case more material is lifted up due 
to the larger energy injection from the photosphere and the 
temperature cannot increase easily owing to the larger ra- 
diative cooling. Hence, the chromosphere extends to higher 
altitude in the active case. Interestingly enough, the obser- 
vation of a young solar-like star, CoRot-2a, by using the 
Rossiter-McLaughlin effect shows that this star possesses 
an extended chromosphere of up to 1.16 times of the stellar 
radius (Czesla et al. 2012). 

Figure[2]displays the fractions of the kinetic energy flux 

2 

of the stellar wind, M^, (left) and the radiative loss from 
the transition region and corona (right) of each run against 
the total injected Poynting energy associated with Alfven 
waves, L wavc , from the photosphere. Both panels show large 
scatters in the vertical direction, reflecting varieties of /. 
The sums of these two fractions are much less than unity in 
all the runs, because most of the input energies are reflected 
back downward. 

The left panel shows that the fraction of the kinetic en- 
ergy flux initially increases with increasing input energy; 
the increase of the kinetic energy of the stellar wind is more 
rapid than linear on the input energy, which is related with 
the global instability owing to the reflection and nonlin- 
ear dissipation of Alfvenwaves as explained above (Suzuki 
2012). However, it eventually saturates and decreases on in- 
creasing iwavc- On the other hand, the radiative loss shows 
saturation as well but the overturning trend is not so distinc- 
tive. 

Combining these two panels, we can get the relation be- 
tween radiative flux and mass loss rate (Fig. |3), which can 
be directly compared with the figures presented in Wood et 
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Fig. 3 Mass loss rate normalized by the present solar 
value on radiative flux from transition region and corona. 

al. (2002, 2005). As expected from Fig. [2] the mass loss rate 
is positively correlated with radiative flux initially, while it 
eventually saturates, or even slightly decreases. The main 
reason of the saturation is that most of the input energies are 
exhausted for the radiative losses besides reflection in these 
saturated cases. Namely, the increase of M is mainly done 
by the increase of density, which enhances radiative losses. 
Finally, no more energy remains for the kinetic energy of 
the stellar wind. We discuss this saturation mechanism in 
more detail in a forthcoming paper (Suzuki et al. 2012, in 
preparation). 

4 Mass loss from red giants 

After the hydrogen in the core is exhausted, a star evolves 
to a red giant stage. Because red giant stars have surface 
convection zones, it is expected that the fluctuations the sur- 
faces would give a significant contribution to the stellar winds. 
On the other hand, the expansion of the radius largely af- 
fects the dynamics through the change of the surface gravity. 
We have carried out MHD simulations for red giant winds 
(Suzuki 2007), of which we briefly summarize the results. 
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Fig. 4 From the top to the bottom, radial outflow velocity, 
v r (km s -1 ), temperature, T (K), and density, p(g cm -3 ), 
are plotted. The dotted, dash-dotted, solid, and dashed lines 
are the results of stellar radii, R = R Q (the present Sun), 
3.1-R© (sub-giant), IORq (red giant), and 31-R© (red giant), 
respectively. 

Figure|4]presents the evolution of stellar winds of a 1M© 
star from main sequence to red giant stages. The surface am- 
plitude is estimated from the convective flux (Renzini et al. 
1977), whereas this should be tested by comparison with re- 
cently observed chromospheric fluxes (e.g. Perez Martinez 
et al. 201 1). The middle panel shows that the average tem- 
perature drops suddenly from T ~ 7 x 10 5 K in the sub- 
giant star (dash-dotted) to T < 10 5 K in the red giant stars, 
which is consistent with the observed "dividing line" (Lin- 
sky & Haisch 1979). The main reason of the disappearance 
of the steady hot coronae is that the sound speed (w 150 
km s _1 ) of ~ 10 6 K plasma exceeds the escape speed, 
Vesc (r) = y/2GAI*/r, at a few stellar radii in the red gi- 
ant stars; the hot corona cannot be confined by the gravity 
any more in the atmospheres of the red giant stars. There- 
fore, the material flows out before heated up to coronal tem- 
perature as an extended chromosphere (Schroder & Cuntz 
2005) with intermittent hot bubbles (Suzuki 2007; Crowley 
et al.2008). In addition, the thermal instability of the radia- 
tive cooling function (Landini & Monsignori-Fossi 1990) 
plays a role in the sudden decrease of temperature; since 
the gas with T = 10 5 — 10 6 K is unstable, the temperature 
quickly decreases from the subgiant to red giants. 

The densities of the winds increase with stellar evolu- 
tion due to the decrease of the surface gravity. The wind 
velocities are much smaller than the surface escape speeds, 



because the onset locations of the winds are around sev- 
eral stellar radii (Harper et al.2009). Through the stellar 
evolution, the mass loss rate increases due to the increase 
of the stellar surface (oc R 2 ) and the increase of the den- 
sity. From the main sequence star to the read giant star with 
i?* = 31i?0, the mass loss rate increases 10 5 — 10 6 times. 
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